Microcanonical Entropy and Dynamical Measure of Temperature for Systems with 

Two First Integrals 
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We consider a generic classical many particle system described by an autonomous Hamiltonian 
H(x 1 , . . . , x N+2 ) which, in addition, has a conserved quantity V(x 1 , . . . , x N+2 ) — v, so that the 
Poisson bracket {H, V} vanishes. We derive in detail the microcanonical expressions for entropy 
and temperature. We show that both of these quantities depend on multidimensional integrals over 
sub-manifolds given by the intersection of the constant energy hyper-surfaces with those defined by 
V(x 1 , . . . , x N + 2 ) = v. We show that temperature and higher order derivatives of entropy are micro- 
canonical observable that, under the hypothesis of ergodicity, can be calculated as time averages of 
suitable functions. We derive the explicit expression of the function that gives the temperature. 
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For an isolated many-body classical system, the ergod- 
icity makes equivalent thermodynamics and dynamics. 
Thus, in these cases, one can measure thermodynamic 
quantities, like temperature and specific heat, as tem- 
poral averages of suitable functions along almost each 
trajectory of a given system, or, equivalently, as an inte- 
gral over its phase space. The opportunity to pass from 
the dynamics to the microcanonical-thermodynamics and 
vice versa, gives the possibility to choose the smarter way 
to measure a given quantity. Very often the calculation 
of thermodynamic quantities in the microcanonical en- 
semble is an impracticable issue, thus one is forced to 
recur to the canonical ensemble, where these measures 
are more easily performed by resorting to numeric simula- 
tions, e.g. by Monte Carlo method. Of course, for a given 
system this step is performable only if there is equiva- 
lence between canonical and microcanonical ensembles in 
the thermodynamic limit. Furthermore, it is worth men- 
tioning that in many cases, the convergence of thermo- 
dynamic quantities to their asymptotic values, is much 
faster if the averages are computed along the dynam- 
ics, rather than by an important sampling of the canonic 
phase-space. For these cases therefore, the dynamics is 
preferred respect to statistics. Furthermore, only for sys- 
tems described by stable and temperated inter-particle 
interaction potentials is guaranteed the equivalence of the 
statistical ensembles in the thermodynamic limit. Nowa- 
days, several of the most intriguing challenge for modern 
science, deals with systems of size intermediate between 
the macroscopic and the microscopic scales. Systems like 
polymers, DNA-helix, proteins, nanosystems, are large 
enough to allow a statistical treatment, but are abso- 
lutely far from the thermodynamic limit. Thus, for these 
systems, ensemble equivalence is hardly verified and one 
has no option but performing dynamical averages. 

It is in this same spirit that Rugh, in [lj , has presented 
a dynamical approach for measuring the temperature of 



a Hamiltonian system in the microcanonical ensemble. 
He has shown that for an ergodic classical system, which 
has only one conserved quantity, i.e. the energy, the in- 
verse temperature 1/T, can be calculated as a tempo- 
ral average of a suitable functional along the dynamics. 
The Rugh's seminal work has stimulated several papers 
which have aroused widespread interest, especiall" 
among those who simulate the properties of liquids p 
The calculation given by Rugh in Ref. [l[ provides a 
microcanonical definition of temperature that allows its 
measure also in systems with nonstandard Hamiltonians. 
Nevertheless, this calculation works for systems with only 
one conserved quantity, i.e. the total energy. In the 
present paper we extend the calculation of entropy and 
of microcanonical temperature to the case of Hamiltonian 
systems with two first integrals of motion. 

In the present paper we consider a classical system of 
(N + 2) degrees of freedom (with TV > 0), which is de- 
scribed by an autonomous Hamiltonian H, and which 
has a further independent conserved quantity V, such 
that {H, V} = 0. We derive in detail the microcanon- 
ical expressions for entropy, i.e. the expression of the 
microcanonical invariant measure, and the temperature, 
moreover we give a formula to derive, recursively, all 
order of derivatives of the entropy. We show that en- 
tropy and temperature depend on multidimensional in- 
tegrals over sub-manifolds given by the intersection of 
the constant energy hyper-surfaces with those defined 
by V{x 1 , . . . , x N+2 ) — v. In particular, we show that 
temperature and higher order derivatives of entropy are 
microcanonical observable that, under the hypothesis of 
ergodicity, can be calculated as time averages of suitable 
functions. In Ref. [f| it has been studied the microcanon- 
ical ensemble of a classical system, whose Hamiltonian is 
parameter dependent and in presence of other first inte- 
grals and it has been showed a method, alternative to the 
present one, that allows one to obtain the first derivative 
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of entropy respect to the conserved quantities. The aim 
of the present paper is to derive explicitly the functional 
by means of which the temperature can be calculated as a 
microcanonical average, in the case of a generic classical 
system, described by a many-body Hamiltonian with one 
further conserved quantity. Furthermore, our method al- 
lowed us to derive an iterative formula that gives the 
derivatives of S(E) of all orders, for this class of systems. 
By this formula, one can measure more general quantities 
like, e.g. the specific heat. 

Let H(x, x N+1 , x N + 2 ) be a classical Hamiltonian de- 
scribing an autonomous many-body system whose co- 
ordinates and canonical momenta (qi , pi , . . . ) are repre- 
sented as (N + 2)-component vectors (x, x N+1 , x N+2 ) € 
K N + 2 , and let V(x, x N+1 , x N + 2 ) be a further conserved 
quantity which is in involution with H. We shall 
assume that V is a smooth function of the coordi- 
nates. The system's motion takes place on the mani- 
folds Ai = fl V u , where the Y, E = {(x,x N+1 ,x N + 2 ) € 
M. N + 2 \H(-x,x N + 1 ,x N+2 ) = E} are energy level sets and 
V u = {(x,x N+1 ,x N + 2 ) e R N+2 \V(x, x N + 1 ,x N + 2 ) = u} are 
subsets of R N + 2 where V is constant. Among the equiv- 
alent expressions allowed for the microcanonical entropy 
S(E), the surface entropy 0] 

S(E) = ln J d N xdx N+1 dx N+2 5{H(x,x N + 1 ,x N+2 )-E)x 

5(V(x,x N+1 ,x N+2 )-u) (1) 

has an interesting and useful geometric interpretation 
that we shall derive, following the calculation shown in 
the chapter on the theory of surfaces of Ref. 0. Con- 
sistently with Ref. Q, we shall assume the level sets of 
H and V to be non-singular hyper-surfaces. Even if the 
energy level surfaces (or the level sets of V) in general 
constitute a singular foliation, thus for some values of E 
(or u) the energy surface (or V u ) is not a differential man- 
ifold, for generic values of E (or u) this is not an issue. 
For a generic point xq € M. of a non-singular level-set 
of H and V, VH(x ) ^ and VV(x ) ^ 0. Further- 
more, since H and V are in involution, almost everywhere 
VH(xq) and W{xq) are independent vectors. Thus, in 
a neighborhood of xq we can suppose of reorder the coor- 
dinate indices in such a way that dH /dx N+1 dV /dx N+2 — 
dV/dx N + 1 dH/dx N + 2 ^ for each x of the neighbor- 
hood. Now, each non-singular manifold M, can be par- 
titioned by a family T of not overlapping subsets [la ]. 
With the further condition that, in each subset we can 
reorder the coordinate indices, as said above, so that 
dH/dx N + 1 dV/dx N + 2 - dV /dx N ^dH /dx N + 2 ^ every- 
where in the subset. Now, by using the first N co- 
ordinates we give a parametric description of the same 
subset. Thus, in each of the subsets of J- we can choose 
f a — id for a — 1, ... , N, and let us set g(y) := (y) 
and h(y) := f N+2 (y), where y E R N . The metric induced 
by R N + 2 on M results 



where <9 Q » = d • /dx a , whereas the its determinant can 
be derived by straightforward calculations and it results 
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V> (d^gd„h - d^hd^gf . (3) 



The derivatives d a g, d a h, can be expressed as follows 

d a g = [d N +2Vd a H - 8n+2 Hd a V] /D (4) 
d a h = [d N +iHd a V - dN+iVd a H] /D , (5) 

where D =d N+1 Hd N+2 V — 8 N+1 Vd N+2 H [16j. From the 
expression above we derive the sub-manifold volume form 



dr 



iff r jJV 
d xJn = d x — 
V / D 
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On the other hand, expression (fTJ) can be cast in the 
following form. For each point x we can introduce the 
following variables change 



y 



1 = h]^{x N +\x N + 2 ) . 



V 



hl(x 



N + l ,j,N + 2 



), (8) 



with the inverse transformations 

x»» =Gl{y\y 2 ), x^ 2 =G 2 x (y\y 2 ), (9) 

where h\{x N + 1 ,x N + 2 ) = H(x, x N+1 , x N + 2 ) - E, 
hl(x N+1 ,x N + 2 ) = V{x,x N +\x N+2 )-u, and G\{y,z) and 
Gx(y, z) are the respective inverse transformations. From 
now on we will suppress the sub- index x . After these def- 
initions Eq. (fTJ) can be expressed as follows 

S(E)=\n J d N xd y 1 dy 2 S(y 1 )6(y 2 )\J(y\y 2 )\, 

where \J{y 1 ,y 2 )\ is the determinant of the Jacobian ma- 
trix d(x N+1 , x N+2 )/d(y 1 , y 2 ). The inverse of the Jacobian 
matrix can be derived by Eqs. ©-(I!]) and it results 



dp,gd v g + d^hd^h , 



(2) 



dxG 1 = -d 2 h 2 /D = -dN + 2V/D 1 
<9iG 2 = d 1 h 2 /D = d N+ iV/D, 
d 2 G 1 = d 2 h}/D = 8n+ 2 H/D, 
d 2 G 2 = -dih 1 /D = -8n + 2H/D . 

Thus, the microcanonical entropy [l7| results 



(10) 

(11) 

(12) 
(13) 



(14) 
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It is worth emphasizing that the invariant measure is in- 
dependent from the partition ofAi, since W is invariant 
under exchange of the indices of the coordinates. 

In order to derive the temperature in the micro- 
canonical ensemble, according to the definition T(E) = 
(dS(E) /dE)^ 1 , we shall use the following generaliza- 
tion [18 of the Federer-Laurence derivation formula 
[l3j]. The flux <I>£ with a non- vanishing component in 
the direction orthogonal to the constant energy hyper- 
surfaces, but tangent to the level hyper-surfaces of V, 
can be defined by the vector £ = n H — (n H ■ n v )n v , 
where n H = \/H/\\ y H\\ and n v = V^VII V ^11- Let us 
define n 5 = £/||£||, so the generalized derivation formula 
results 



Qk 

dE* 



drip 



M 



drA k (VO 



(15) 



M 



where [lj| 
A(.) 



1 



{n v • V)(« 5 )] • (16) 



The proof of the extention of the Federer-Laurence theo- 
rem to varieties of codimension two, is rather complicated 
and lengthy. Furthermore it is outside of the main mo- 
tivation of the present paper, thus it will be given in a 
further paper [14[. After Eqs. (fTT|) and (IT51) . we obtain 
that the inverse temperature is given by 



and is defined for D ^ (z ^ w) that is on the subsets 



{{x,y,±\l ^ >T 



1 — x 2 — y 2 



y 2 < !}• 



The sum of the integrals of y/rj upon these subsets brings 
to the result 4-7T which is the right hyper-surface vol- 
ume. This solution can be checked noting that in the 
original problem the condition V = z + w = can be 
lifted by a suitable change of variables. Let us intro- 
duce the variables s — (z — w )/V2 and t = (z + w)/y/2. 
H and V can be expressed in the new variables as 
H = x 2 + y 2 + s 2 + 1 2 = 1 and V = t = 0. Thus the met- 
ric r] is that one of the three-dimensional unitary sphere 
x 2 + y 2 + s 2 = 1 that is (see Q) r] = 1 + (x/s) 2 + (y/s) 2 . 



This indeed is the expression (|20|) . The sum of the inte- 
grals of yjr\ upon the two hemispheres, s > and s < 0, 
gives the right value Air. 

ii) Let us now consider a simple case where all the ge- 
ometric quantities can be explicitly calculated in order 
to check the Eqs. (fT5")) and (fTB]) . We shall consider 



H 



X o + X ?. 



E, a sphere of radius \f~E~ in three 
dimensions and V — x^/ \/ x\ + x\ = u, that is a cone 
with an angle arctan(u) at the vertex. In this case Ai is 
a circle of radius a = y/E/ (1 + u 2 ). 

If we choose ip = 1 , we get easly 



1 



T(E,V) 



where 



w 



X/H ■ n« 



= <*(*)) 



(n v ■ y)(rce) 
W 



(17) 



(18) 



and () stands for the microcanonical average. 

When (M. N+2 ,H) is ergodic with respect to the Liou- 
villc measure dr/W restricted to a nonsingular manifold 
M. for almost every initial condition x(0) € Ai one has 
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M 



d_ 

dE' 



2na 



~E 



The terms that appear in Eq. (JTHJ) result sjH-n^ = 2y/E, 

V" 5 = 2/y/E and n v ■ (n v ■ v)( n? ) = 1/Ve~- Thus 
A(l) = 1/(2E) and consequently J M drA (1) = ira/E. 



By choosing ip 
dr 



d_ 
~dE 



1/W we get W = 2E/a 2 , and 
d a 2 



M 



W 3E2E 



-2ira 



2(l + u 2 ) 3 / 2 VE 



(21) 



T(E,V) 



lim — 

s— >oo S 



ds' [$(«(«'))] 



(19) 



Simple geometric applications, 
i) As a simple application we shall derive the invariant 
metric n in the case of a simple form for H as that one 
of a four-dimensional hypersphere H = x 2 + y 2 + z 2 + w 2 
of unit radius, and with a condition given by the hyper- 
plain V = z + w = 0. By direct calculations we derive 
d x H = 2x, d y H = 2y, d z H = 2z, d w H = 2w, thus it 
results D = d z Hd w V - d z Vd w H = 2(z-w). With the 
notations introduced above, Eqs. ((4]) and (J5)) become 
9x9 = x/(z - w), d y g = y/(z- w), d x h = -x/(z - w), 
-y I ' {z — w) , thus the invariant metric results 



and d y h 



II 



(20) 



Again, Eq. ([H]) contains \/H ■ n« = 2\[E, \j{n^ /W) = 
a 2 /E 3 / 2 and n v ■ (n v ■ Sj){m)/W = a 2 /(2E 3 ^ 2 ). Thus 
A(l/W) = a 2 /(4E 2 ) and consequently J M drA (1/W) = 
ira 3 /(2E 2 ) which indeed coincides with result (l2~Tj) . 

Dynamical system. 
Let us now consider a lattice system described by the 
Hamiltonian 



H 



,2 \2 



(p m p m +i + q m q m +i) (22) 



and the usual Poisson brackets {g m ,Pn} — 5 m .n for 
n,m — 1,...,M, with periodic boundary conditions. 
The dynamics generated by this Hamiltonian conserves 
the quantity N = E TO («m +Pm)A thu s the E qs- CLH) 
and (|18p are in order to calculate the microcanonical tem- 
perature. The ground-state is achieved by solving the 
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equation 5(H — \N) = in which the Lagrangian multi- 
plier x has been introduced to take in account the conser- 
vation of N. By direct calculations, we got the solution 
q 0rn = y/2N/M := q , Po m = and X = vN/M - 2 
20]. Small fluctuations around this ground-state cor- 
respond to T > 0, in the following we show that this 
prevision is verified by Eq. (TIT)) and (TT51) . whereas the 
formula for 1/T given in that holds in the case of 
systems with only one first integral (energy), does not 
work. The reason of this comparison is to show that 
the equation derived in [lj, which is valid in the case 
of system with only one first integral, cannot be used 
in the case of systems with more than one first inte- 
grals In Ref. p has been derived an equation sim- 
ilar to Eqs. (TIT)) and (fT5)l . The inverse temperature 
1/T can be derived by Eqs. (12) and (17) of Ref. @ 
after having found two vectors, Xq and X\, such that 
dH(Xo) = 1, dH{X x ) = 0, dV(X ) = and dV{Xx) = 
1. One can use Xo — cvv^H/d — cyjyVV/d and 
Xi = -c VH VH/d + c HH VV/d, where c HH = ||ViJ|| 2 , 
cvv = || W|| 2 , c VH = WViJ, and d = c HH c V v ~ c 2 VH . 
Thus it results 1/T = (V • Xq) , The latter expression 
seems inequivalent to Eqs. (fl~7|) (fT8|) . Indeed, it does not 
contain 1/W which is related to the invariant measure. 
It is worth emphasizing that an analogous term l/||iJ||, 
which is related to the microcanonical measure of a sys- 



tem with one first integral, indeed appears in the formula 
for I jT derived in [l| . In any case a comparison between 
these two equations, would require a numerical simula- 
tion, but this is out of the aim of the present paper. By 
expanding H in terms of the displacements Q n and P n 
from the minimum points q n = qg and p n — 0, and by cal- 
culating the terms appearing in (|I7[) and (|I8[) in the limit 
Qn, Pn — > 0, after some boring algebra we got the correct 
low-energy temperature, being 1/T = ( C E > ) MC —¥ oo. 

It is worth emphasizing that by using the formula de- 
rived in Ref. [J 1/T = (V (Vi?/|| VH || 2 ) ) which is 
correct for systems with only one first integral, one ob- 
tains the erroneous result: T = [(f 9 2 — 2) q$\ /(2v(^). 

In conclusion, we have presented a dynamical approach 
for measuring the temperature of a Hamiltonian system 
with two first integrals in the microcanonical ensemble. 
The formula we have derived allows one to measure the 
inverse temperature as a time-average, instead of as an 
average over the phase-space, also in the case of systems 
with two first integrals. Furthermore, by Eq. i U5)) higher 
orders of derivatives of S(E) can be obtained. Therefore 
the method here presented allows one to measure g.e. the 
specific-heat. 

I'm indebted with Prof. M. Pettini, Prof. G. Vezzosi 
and Prof. R. Livi for the usefull discussions. 
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